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(57) Abstract 

A technique is described for identifying mutations, if 
any, present in a biological sample, from a pre-selected 
set of known mutations. The method can be applied to 
DNA, RNA and peptide nucleic acid (PNA) microarrays. 
The method analyzes a dot spectrogram representative of 
quantized hybridization activity of oligonucleotides in the 
sample to identify the mutations. In accordance with 
the method, a resonance pattern is generated which is 
representative of nonlinear resonances between a stimulus 
pattern associated with the set of known mutations and the 
dot spectrogram. The resonance pattern is interpreted to a 
yield a set of confirmed mutations by comparing resonances 
found therein with predetermined resonances expected for 
the selected set of mutations. In a particular example, the 
resonance pattern is generated by iteratively processing the 
dot spectrogram by performing a convergent reverberation 
to yield a resonance pattern representative of resonances 
between a predetermined set of selected Quantum Express or 
Functions and the dot spectrogram until a predetermined 
degree of convergence is achieved between the resonances 
found in the resonance pattern and resonances expected for 
the set of mutations. The resonance pattern is analyzed 
to a yield a set of confirmed mutations by mapping the 
confirmed mutations to known diseases associated with the 
p re-selected set of known mutations to identify diseases, 
if any, indicated by the biological sample. By exploiting 
a resonant interaction, mutation signatures may be robustly 
identified even in circumstances involving low signal to noise ratios or, in some cases, negative signal to noise ratios. 



Blxat Output from 
CU»»hc*l Proc*»m 





FOR THE PURPOSES OF INFORMATION ONLY 
Codes used to identify States party to the PCT on the front pages of pamphlets publishing international applications under the PCT. 



AL 


Albania 


ES 


Spain 


LS 


Lesotho 


SI 


Slovenia 


AM 


Armenia 


FI 


Finland 


LT 


Lithuania 


SK 


Slovakia 


AT 


Austria 


FR 


France 


LU 


Luxembourg 


SN 


Senegal 


AU 


Australia 


GA 


Gabon 


LV 


Latvia 


SZ 


Swaziland 


AZ 


Azerbaijan 


GB 


United Kingdom 


MC 


Monaco 


TD 


Chad 


BA 


Bosnia and Herzegovina 


GE 


Georgia 


MD 


Republic of Moldova 


TG 


Togo 


BB 


Barbados 


GH 


Ghana 


MG 


Madagascar 


TJ 


Tajikistan 


BR 


Belgium 


GN 


Guinea 


MK 


The former Yugoslav 


TM 


Turkmenistan 


BF 


Burkina Faso 


GR 


Greece 




Republic of Macedonia 


TR 


Turkey 


BG 


Bulgaria 


HU 


Hungary 


ML 


Malt 


TT 


Trinidad and Tobago 


BJ 


Benin 


IE 


Ireland 


MN 


Mongolia 


UA 


Ukraine 


BR 


Brazil 


IL 


Israel 


MR 


Mauritania 


UG 


Uganda 


BY 


Belarus 


IS 


Iceland 


MW 


Malawi 


US 


United States of America 


CA 


Canada 


IT 


Italy 


MX 


Mexico 


UZ 


Uzbekistan 


CF 


Central African Republic 


JP 


Japan 


NE 


Niger 


VN 


Viet Nam 


CG 


Congo 


KE 


Kenya 


NL 


Netherlands 


YU 


Yugoslavia 


CH 


Switzerland 


KG 


Kyrgyzstan 


NO 


Norway 


ZVV 


Zimbabwe 


CI 


Cdte d'lvoire 


KP 


Democratic People's 


NZ 


New Zealand 






CM 


Cameroon 




Republic of Korea 


PL 


Poland 






CN 


China 


KR 


Republic of Korea 


PT 


Portugal 






CU 


Cuba 


KZ 


Kazakstan 


RO 


Romania 






cz 


Czech Republic 


LC 


Saint Lucia 


RU 


Russian Federation 






DE 


Germany 


U 


Liechtenstein 


SD 


Sudan 






OK 


Denmark 


LK 


Sri Lanka 


SE 


Sweden 






EE 


Estonia 


LR 


Liberia 


SG 


Singapore 







WO 00/52625 PCTAJSOO/04076 
METHOD AND APPARATUS FOR ANALYZING HYBRIDIZED BIOCHIP 
PATTERNS USING RESONANCE INTERACTIONS EMPLOYING QUANTUM 

EXPRESSOR FUNCTIONS 

5 Field of the Invention 

The invention generally relates to techniques for analyzing biological samples such as 
DNA, RNA, or protein samples and in particular to techniques for analyzing the output patterns of 
hybridized biochip microarrays. 

Background of the Invention 

10 A variety of techniques have been developed to analyze DNA or other biological samples 

to identify diseases, mutations, or other conditions present within a patient providing the sample. 
Such techniques may determine, for example, whether the patient has any particular disease such 
as cancer or AIDS, or has a predisposition toward the disease, or other medical conditions present 
in the patient.DNA-based analysis may be used either as an in-vitro or as an in-vivo control 

1 5 mechanism to monitor progression of disease, assess effectiveness of therapy or be used to design 
dosage formulations. DNA-based analysis is used verify the presence or absence of expressed 
genes and polymorphisms. 

One particularly promising technique for analyzing biological samples uses a DNA-based 
microarray (or microelectronics biochip) which generates a hybridization pattern representative of 

20 the characteristics of the DNA within the sample. Briefly, a DNA microarray includes a 
rectangular array of immobilized single stranded DNA fragments. Each element within the array 
includes few tens to millions of copies of identical single stranded strips of DNA containing 
specific sequences of nucleotide bases. Identical or different fragments of DNA may be provided 
at each different element of the array. In other words, location (1,1) contains a different single 

25 stranded fragment of DNA than location (1,2) which also differs from location (1,3) etc. Certain 
biochip designs may replicate the nucleotide sequence in multiple cells. 

DNA-based microarrays deploy chemiluminiscence, fluorescence or electrical 
phenomenology to achieve the analysis. In methods that exploit fluorescence imaging, a target 
DNA sample to be analyzed is first separated into individual single stranded sequences and 

30 fragmented. Each sequence being tagged with a fluorescent marker molecule. The fragments are 
applied to the microarray where each fragment binds only with complementary DNA fragments 
already embedded on the microarray. Fragments which do not match any of the elements of the 
microarray simply do not bind at any of the sites of the microarray and are discarded during 
subsequent fluidic reactions. Thus, only those microarray locations containing fragments that bind 

1 
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complementary sequences within the target DNA sample will receive the fluorescent molecules. 
Typically, a fluorescent light source is then applied to the microarray to generate a fluorescent 
image identifying which elements of the microarray bind to the patient DNA sample and which do 
not. The image is then analyzed to determine which specific DNA fragments were contained 
5 within the original sample and to determine therefrom whether particular diseases, mutations or 
other conditions are present in the patient sample. 

For example, a particular element of the microarray may be exposed to fragments of DNA 
representative of a particular type of cancer. If that element of the array fluoresces under 
fluorescent illumination, then the DNA of the sample contains the DNA sequence representative of 
10 that particular type of cancer. Hence, a conclusion can be drawn that the patient providing the 
sample either already has that particular type of cancer or is perhaps predisposed towards that 
cancer. As can be appreciated, by providing a wide variety of known DNA fragments on the 
microarray, the resulting fluorescent image can be analyzed to identify a wide range of conditions. 
Unfortunately, under conventional techniques, the step of analyzing the fluorescent pattern 
1 5 to determine the nature of any conditions characterized by the DNA is expensive, time consuming, 
and somewhat unreliable for all but a few particular conditions or diseases. One major problem 
with many conventional techniques is that the techniques have poor repeatability. Hence, if the 
same sample is analyzed twice using two different chips, different results are often obtained. Also, 
the results may vary from lab to lab. Consistent results are achieved only when the target sample 
20 has high concentrations of oligonucleotides of interest. Also, skilled technicians are required to 
prepare DNA samples, implement the hybridization protocol, and analyze the DNA microarray 
output possibly resulting in high costs. One reason that repeatability is poor is that the signatures 
within the digitized hybridization pattern (also known as a "dot spectrogram") that are 
representative of mutations of interest are typically very weak and are immersed in considerable 
25 noise. Conventional techniques are not particularly effective in extracting mutation signatures from 
dot spectrograms in low signal to noise circumstances. Circumstances wherein the signal to noise 
ratio is 0 to strongly negative (-2 to -30dB) are particularly intractable. 

Accordingly, it would highly desirable to provide an improved method and apparatus for 
analyzing the output of the DNA microarray to more expediently, reliably, and inexpensively 
30 determine the presence of any medical conditions or concerns within the patient providing the 
DNA sample. It is particularly desirable to provide a technique that can identify mutation 
signatures within dot spectrograms even in circumstance wherein the signal to noise ration is 
extremely low. It is to these ends that aspects of the invention are generally drawn. 

Referring now to FIG. 1, conventional techniques for designing DNA microarray chips and 

2 
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for analyzing the output thereof will now be described in greater detail. Initially, at step 100, 
fluorescently labeled primers are prepared for flanking loci of genes of interest within the DNA 
sample. The primers are applied to the DNA sample such that the fluorescently labeled primers 
flank genes of interest. At step 102, the DNA sample is fragmented at the locations where the 
5 fluorescently labeled primers are attached to the genes of interest to thereby produce a set of DNA 
fragments, also called "oligonucleotides" for applying to the DNA microarray. 

In general, there are two types of DNA microarrays: passive hybridization microarray s and 
active hybridization microarTays. Under passive hybridization, oligonucleotides characterizing the 
DNA sample are simply applied to the DNA microarray where they passively attach to 

1 0 complementary DNA fragments embedded on the array. With active hybridization, the DNA array 
is configured to externally enhance the interaction between the fragments of the DNA samples and 
the fragments embedded on the microarray using, for example, electronic techniques. Within FIG. 
1, both passive hybridization and active hybridization steps are illustrated in parallel. It should be 
understood that, currently for any particular microarray, either the passive hybridization or the 

15 active hybridization steps, but not both, are typically employed. Referring first to passive 
hybridization, at step 104 a DNA microarray chip is prefabricated with oligonucleotides of interest 
embedded or otherwise attached to particular elements within the microarray. At step 106, the 
oligonucleotides of the DNA sample generated at step 102 are applied to the microarray. 
Oligonucleotides within the sample that match any of the oligonucleotides embedded on the 

20 microarray passively bind with the oligonucleotides of the array while retaining their fluorescently 
labeled primers such that only those locations in the microarray having corresponding 
oligonucleotides within the sample receive the primers. It should be noted that each individual 
nucleotide base within the oligonucleotide sequence (with lengths ranging from 5 to 25 base pairs) 
can bond with up to four different nucleotides within the microarray, but only one oligonucleotide 

25 represents an exact match. When illuminated with fluorescent light, the exact matches fluoresces 
most effectively and the non-exact matches fluoresce considerably less or not at all. 

At step 108, the DNA microarray with the sample loaded thereon is placed within a fluidics 
station provided with chemicals to facilitate the hybridization reaction, i.e., the chemicals facilitate 
the bonding of the oligonucleotide sample with corresponding oligonucleotides within the 

30 microarray. At step 110, the microarray is illuminated under fluorescent light, perhaps generated 
using an ion-argon laser, and the resulting fluorescent pattern is digitized and recorded. 
Alternately, a photograph of the fluorescent pattern may be taken, developed, then scanned into a 
computer to provide a digital representation of the fluorescent pattern. In any case, at step 112, the 
digitized pattern is processed using dedicated software programs which operate to focus the digital 
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pattern and to subsequently quantize the pattern to yield a fluorescent intensity value for each array 
within the microarray pattern. At step 114, the resulting focused array pattern is processed using 
additional software programs which compute an average intensity value at each array location and 
provides for necessary normalization, color compensation and scaling. Hence, following step 1 14, 
5 a digitized fluorescent pattern has been produced identifying locations within the microarray 
wherein oligonucleotides from the DNA sample have bonded. This fluorescent pattern is referred 
to herein as a "dot spectrogram". 

In existing biochips that actively initiate, facilitate or selectively block hybridization, a 
DNA microarray is prefabricated for active hybridization at step 116. At step 118, the DNA 

10 sample is applied to the active array and electronic signals are transmitted into the array to help 
facilitate bonding between the oligonucleotides of the sample and the oligonucleotides of the array. 
The microarray is then placed within a fluidics station which further facilitates the bonding. 
Thereafter, at step 122, an electronic or fluorescent readout is generated from the microarray. 
When electrical output signals from the biochip array are used to quantify and classify the post- 

15 hybridization output, the output signal is indicative of the number oligonucleotide fragments 
bonded to each site within the array. At step 124 the electronic output is processed to generate a 
dot spectrogram similar or identical to the dot spectrogram generated using the optical readout 
technique of steps 110-114. Hence, regardless of whether steps 104-114 are performed or steps 
116-124 are performed the result is a dot spectrogram representative of oligonucleotides present 

20 within the DNA sample. Here it should be noted that some conventional passive hybridization 
DNA microarrays provide electronic output and some active hybridization microelectronic arrays 
provide optical readout. Thus, for at least some techniques, the output of step 108 is processed in 
accordance with steps 122 and 124. For other techniques, the output of step 120 is processed in 
accordance with steps 110-114. Again, the final results are substantially the same, i.e., a dot 

25 spectrogram. 

At step 126, the dot spectrogram is analyzed using clustering software to generate a gene 
array amplitude readout pattern representative of mutations of interest within the target DNA 
sample. In .essence, step 126 operates to correlate oligonucleotides represented by the dot 
spectrogram with corresponding DNA mutations. Next, at step 128, the resulting digital 
30 representation of the mutations of interest are processed using mapping software to determine 
whether the mutations are representative of particular diagnostic conditions, such as certain 
diseases or conditions. Hence, step 128 operates to perform a mutation-to-diagnostic analyses. 
Finally, at step 130 the diagnostic conditions detected using step 128 are evaluated to further 
determine whether or not the diagnostic, if any, can properly be based upon the DNA sample. 
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Classical methods such as probabilistic estimator such as minimum a posteriori (MAP) estimator, 
maximum likelihood estimator (MLE) or inferencing mechanism may be used to render a 
diagnostic assessment. 

As noted above, it would be desirable to provide improved techniques for analyzing the 
5 outputs for DNA microarrays to more quickly, reliably and inexpensively yield a valid diagnostic 
assessment To this end, the invention is directed primarily to providing a sequence of steps for 
replacing steps 114-130 of FIG. 1. 

Summary of the Invention 

10 In accordance with a first aspect of the invention, a method is provided for analyzing an 

output pattern of a biochip to identify mutations, if any, present in a biological sample applied to 
the biochip. In accordance with the method, a resonance pattern is generated which is 
representative of resonances between a stimulus pattern associated with a set of known mutations 
and the output pattern of the biochip. The resonance pattern is interpreted to yield a set of 

15 confirmed mutations by comparing resonances found therein with predetermined resonances 
expected for the selected set of mutations. 

In an exemplary embodiment, the biological sample is a DNA sample and the output 
pattern being analyzed is a quantized dot spectrogram generated by a hybridized oligonucleotide 
microarray. The resonance pattern is generated by iteratively processing the dot spectrogram by 

20 performing a convergent reverberation to yield a resonance pattern representative of resonances 
between a predetermined set of selected Quantum Expressor Functions and the dot spectrogram 
until a predetermined degree of convergence is achieved between the resonances found in the 
resonance pattern and resonances expected for the set of mutations. The resonance pattern is 
analyzed to yield a set of confirmed mutations by mapping the confirmed mutations to known 

25 diseases or diagnostic conditions of interest, associated with the pre-selected set of known 
mutations to identify diseases, if any, indicated by the DNA sample. A diagnostic confirmation is 
then made by taking the identified diseases and solving in reverse for the associated Quantum 
Expressor Functions and then comparing those Quantum Expressor Functions with ones expected 
for the mutations associated with the identified disease to verify correspondence. If no 

30 correspondence is found, a new sub-set of known mutations are selected and the steps are repeated 
to determine whether any of the new set of mutations are present in the sample. 

In the exemplary embodiment the set of nonlinear Quantum Expressor Functions are 
generated are follows. A set of mutation signatures representative of the pre-selected set of known 
mutations is input A representation of a microarray oligonucleotide pattern layout for the 
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microarray, from which the dot spectrogram is generated, is also input. Then a set of resonant 
interaction parameters are generated which are representative of mutation pattern interactions 
between elements of the microarray including interactions from a group including element-to- 
element interactions, element-to-ensemble interactions, ensemble-to-element interactions, and 
5 ensemble-to-ensemble interactions. Then the set of nonlinear Quantum Expressor Functions are 
generated from the set of resonant interaction patterns by matching selected harmonics of the 
power spectral density (PSD) amplitude of a coded mutation signature, corresponding to the pre- 
selected mutation set of interest, to that of a pre-determined quantum-mechanical Hamiltonian 
system so that stochastic and deterministic time scales match, and the time scales couple back to 

1 0 noise statistics and degree of asymmetry. 

Also in the exemplary embodiment, the dot spectrogram is differentially enhanced prior to 
the generation of the resonance pattern by refocusing the dot spectrogram to yield a re-focused dot 
spectrogram; cross-correlating the re-focused dot spectrogram; applying a local maxima filter to 
the correlated re-focused dot spectrogram to yield a maximized dot spectrogram; re-scaling the 

1 5 maximized dot spectrogram to yield a uniformly re-scaled dot spectrogram; and then re-scaling the 
uniformly re-scaled dot spectrogram to amplifying local boundaries therein to yield a globally re- 
scaled dot spectrogram. 

By exploiting a resonant interaction, mutation signatures may be identified within a dot 
spectrogram even in circumstances involving low signal to noise ratios or, in some cases, negative 

20 signal to noise ratios. By permitting the mutation signatures to be identified in such 
circumstances, the reliability of dot spectrogram analysis is thereby greatly enhanced. With an 
increase in reliability, costs associated with performing the analysis are decreased, in part, because 
there is less of a requirement for skilled technicians. Other advantages of the invention arise as 
well. 

25 In accordance with a second aspect of the invention, a method of generating a set of 

nonlinear Quantum Expressor Functions is provided. The method includes the steps of inputting a 
set of mutation signatures representative of the pre-selected set of known mutations and inputting a 
representation of a biochip layout. The method also includes the steps of generating a set of 
resonant interaction parameters representative of mutation pattern interactions between elements of 

30 the microarray including interactions from a group including element-to-element interactions, 
element-to-ensemble interactions, ensemble-to-element interactions, and ensemble-to-ensemble 
interactions and generating the set of nonlinear Quantum Expressor Functions from the set of 
resonant interaction patterns. 

Among other applications, principles of the invention are applicable to the analysis of 

6 
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various arrayed biomolecular, ionic, bioelectronic, biochemical, optoelectronic, radio frequency 
(RF) and electronic microdevices. Principles of the invention are particularly applicable to 
mutation expression analysis at ultra-low concentrations using ultra-high density passive and/or 
active hybridization DNA-based microarrays. Techniques implemented in accordance with the 
5 invention are generally independent of the physical method employed to accumulate initial 
amplitude information from the bio-chip array, such as fluorescence labeling, charge clustering, 
phase shift integration and tracer imaging. Also, principles of the invention are applicable to 
optical, optoelectronic, and electronic readout of hybridization amplitude patterns. Furthermore, 
principles of the invention are applicable to molecular expression analysis at all levels of 
1 0 abstraction: namely DNA expression analysis, RNA expression analysis, protein interactions and 
protein - DNA interactions for medical diagnosis at the molecular level. 
Apparatus embodiments are also provided. 

Brief Description of the Drawings 
The features, objects, and advantages of the present invention will become more apparent 
from the detailed description set forth below when taken in conjunction with the drawings in which 
like reference characters identify correspondingly throughout and wherein: 

Fig. 1 is a flow chart illustrating conventional passive and active hybridization DNA 
microarray analysis techniques. 

Fig. 2 is a flow chart illustrating an exemplary method for analyzing the output of a 
hybridized DNA microarray in accordance with the invention. 
Fig. 3 graphically illustrates the method of Fig. 2. 

Fig. 4 is a flow chart illustrating an exemplary method for generating Quantum Expressor 
Functions for use with the method of Fig. 2. 

Fig. 5 is a flow chart illustrating an exemplary method for preconditioning the output of a 
hybridized DNA microarray for use with the method of Fig. 2. 

Detailed Description of Exemplary Embodiments 

With reference to the remaining figures, exemplary embodiments of the invention will now 
30 be described. The invention will be described primarily with respect to an exemplary method for 
analyzing mutations signatures within output patterns of hybridized microarrays generated using 
DNA samples, but principles of the invention apply to the analysis of other protein-based samples 
or to other types of output patterns as well. 
Overview 

7 
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Briefly, the exemplary method exploits, among other features: (a) a novel representation, 
interpretation and mathematical model for the immobilized oligonucleotide hybridization patterns, 
represented via a dot spectrogram; (b) a new "active" biomolecular target detection and 
discrimination method based on quantum resonance interferometry, and (c) a new spatial hashing 
5 function that yields accurate diagnostic assessment 

To this end, the exemplary method exploits a fundamentally different computational 
paradigm for mutation expression detection in pre-enhanced dot spectrogram realizations. The 
method is an innovative modification to dynamically arrayed quantum stochastic resonance (QSR) 
for discrete system analysis. The arraying strategy is a function of the expression pathway of 

10 interest. The method depends on the molecular diagnostic spectrum being addressed. Coupled 
quantum resonators are employed to significantly enhance signal-to-noise (SNR) performance and 
fuse multiple synthetic renormalized dot spectrogram realizations to better detect prespecified 
biomolecular expression patterns. 

Moreover the exemplary method exploits an enhancement in previous extensions to 

1 5 classical stochastic resonance (SR) and array enhanced SR (AESR) in signal processing and sensor 
data analysis. Stochastic resonance is a phenomenon wherein the response to a sensor, modeled in 
terms of a bistable nonlinear dynamical system, is enhanced by applying a random noise element 
and a periodic sinusoidal forcing function. SR occurs when the SNR passes through a maximum as 
the noise level is increased. 

20 Thus an important aspect of the exemplary method involves the coupling of transformed 

and preconditioned discrete microarray outputs to a mathematical model for a quantum-mechanical 
dynamical system with specific properties. When driven in a particular manner, the coupled system 
exhibits a nonlinear response that corresponds to detection of phenomena of interest. The method 
exploits modulation of observables from a "base" (canonical continuous dynamical system), so 

25 that a selected set of spectral properties match a similar selected spectral properties of a discrete 
spatial tessellation substructure from an amplitude spectrogram derived from bioelectronic 
observables. The method further exploits the concept of convolving a discrete spatial system 
(derived from base mutants of interest) with a continuous asymmetric temporal system to derive a 
spatiotemporal input to further convolve with another discrete spatial projection (of an inherently 

30 partially stabilized spatiotemporal system). 

Hence key components of the exemplary biomolecular detection method are: (i) selection 
of a basis system; (ii) generation of designer Quantum Expressor Function (QEF) for coupling 
with the substrate to be analyzed; (iii) generation of a Hamiltonian to describe relaxation dynamics 
of the coupled system; (i v) modulation of resonance parameters to enforce early resonance; (v) and 
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exploitation of resonance suppressors to verify detection. 
The Method 

Referring first to FIG, 2, initially at step 200, a set of mutations of interest is selected. The 
mutations, for example, may be mutations relevant to cancer, AIDS, or other diseases or medical 
5 conditions. At step 201, preconditioner transforms are generated based upon the selected set of 
mutations. The preconditioner transforms are provided to convert mutation nucleotide sequences 
into expected amplitude patterns in the prespecified microarray representation, given a particular 
biochip layout. At step 202, Quantum Expressor Functions are generated based upon the 
Hamiltonian of a pre-selected basis system. The Quantum Expressor Functions are designed to 

10 couple the Hamiltonian for the selected basis system to a predetermined DNA microarray 
configuration to permit a resonance interaction involving the output of the DNA microarray. 
Resonance stimulus is generated, at step 204, using the Quantum Expressor functions. 

What has been summarized thus far are preliminary steps performed off-line for setting up 
the Quantum Expressor Functions and the corresponding resonance stimulus. These steps need be 

1 5 performed only once for a given set of mutations and for a given DNA microarray configuration. 
Thereafter, any number of output patterns from the DNA microarray may be processed using the 
Quantum Expressor Functions to identify whether any of the mutations of the pre-selected set of 
mutations are found therein. Preferably, Quantum Expressor Functions are pre-generated for a 
large set of mutations and for a large set of DNA microarray patterns such that, for each new DNA 

20 microarray output pattern from each new patient sample, the presence of any of the mutations can 
be quickly identified using the predetermined set of Quantum Expressor Functions. In general, the 
aforementioned steps need be repeated only to update the Quantum Expressor Functions to 
accommodate new and different DNA microarray patterns or to if new mutations of interest need 
to be considered. 

25 At step 206, an output pattern (referred to herein as a Dot Spectrogram) is generated using 

a DNA microarray for which Quantum Expressor Functions have already been generated. At step 
208, the dot spectrogram is preconditioned to yield a dot spectrogram tesselation (DST) to permit 
exploitation of a resonance interaction between the dot spectrogram and the Quantum Expressor 
Functions. The actual resonant interaction, which involves convergent reverberations, is 

30 performed at step 210 until a pre-determined degree of convergence is achieved. Once 
convergence is achieved, a resulting resonance pattern is processed at step 212 to identify any 
mutations represented thereby. As will be described below, step 212 is rendered trivial by virtue of 
the aforementioned resonant interaction which is based upon Quantum Expressor Function already 
correlated with the pre-selected mutations. Hence, no complicated analysis is required to interpret 
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the resonance pattern to identify the mutations. Next, at step 214, the mutations are mapped to 
corresponding diseases and conditions to thereby identify any diseases or conditions that the 
patient providing the sample being analyzed is afflicted with. Again, this is a fairly trivial step. 
Finally, at step 216, diagnostic confirmation is preformed to verify that the diseases or conditions 
5 are present in the sample. This is achieved by starting with the found diseases or conditions and 
then performing the steps of the method in reverse. 

FIG. 3 graphically illustrates the operation of the method of FIG. 2 whereby a Quantum 
Expressor Function 300 is generated based on a mutation set 302. A dot spectrogram 304 is 
applied to the Quantum Expressor Function via an interferometric resonance interaction 306 
1 0 yielding a resonance pattern (not shown) from which mutations signatures 308 (representative of 
mutations present in the sample from which the dot spectrogram was generated) may be identified 
by comparison with mutation set 302. 

Now details of the various steps of the method will be provided. 

Mutations Sets 

1 5 Referring again to FIG. 2, the mutation set of interest generated at step 200 is selected by 

identifying oligonucleotides representative of the {Z} mutations of interest. Each oligonucleotide 
is represented by j) which is given by [otocc,. . .a k ], where a = { A,C,T,G} base associated with 
each array cell [a,b] where 10 < k < 25 . The entire set of unique oligonucleotides denoting 
mutations of interest, A(l,m), is given by [5 0 5, 5 k ] where 5 = {A,C,T,G} length | A | = length | \y | , 

20 and 0<|| A-\|/||<k, and the designed in v|/(l,m) oligonucleotide sequence is a perfect complement 
to only A(l,m) for all l,m. 

As part of step 200, an oligonucleotide table is generated which contains the 
oligonucleotide sequences associated with each mutation of interest identified by row and column 
location (ij). The oligonucleotide table is provided for subsequent use at step 212 to map 

25 locations within the dot spectrogram wherein resonance occurs at step 21 0 to oligonucleotides such 
that mutations present in a sample being analyzed are easily identified. Also as part of step 200, a 
mutation table is generated which contains the diseases associated with each mutation of interest. 
The mutation table is provided for subsequent use at step 214 to map mutations identified at step 
212 to specific diseases or other medical conditions such that the diseases can be easily identified. 

30 Preconditioner Transforms 

Given the microarray layout and the mutations of interest (in the coordinate system 
specified by the prespecified exemplary microarray layout) preconditioner transforms are 
employed at step 201 to generate an amplitude pattern and associated point spread function at a 
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signal-to-noise ratio (SNR)=0. 

In the exemplary embodiment, the DNA microarray is an N by M DNA chip array wherein 
an element of the array is referred to herein as an "oxel": o(iJ). 

The pre-hybridization microarray (PEBC) is expressed as: 

N AS 

5 PEBC = X X °&f) • where N and M refer to the linear (row and column) 

i i 

dimensions of the 2-D microarray. 

An inverse Dirichlet Tessellation (IDT) on the PEBC is applied such that a singular value 
decomposition (SVD) on the resulting array yields the location of the mutation of interest at the 
SNR=0 condition. Conceptually the process corresponds to setting up a concave diffusion front at 
1 0 the mutation-centered oxel. 

The equations used are 

15 a f =k a (Dtf*/* where a c = 4/7. 

N f = k N L(Dtf% where a N = 3/7. 
D H = 7/4. 

20 The constants k,, ko and k N are respectively set to 0.856, 0.68, and 1 .34, and for biochip 

dimensions N,M>100, t is typically set to 10 2 * a 2 /D where a = inter-oxel distance and D is a 
coefficient for percolation of the point spread. D is chosen as 0.001 . These equations impose a 
specific IDT, i.e., rapid diffusion on the mutation-centered oxel, and synthetically tessellate the 
biochip. 

25 

Basis System for Quantum Expressor Functions 

To generate Quantum Expressor Functions (QEF) at step 202, a basis system for the QEF is 
first selected. The selection of the basis system and the generation of the QEF's based thereon 
depends, in part, and the characteristics of the DNA microarray. 
30 The numeric value associated with each oxel is given by: 

6{iJ) = a k A k ~* +a, H .4*- 2 + ...+a 1 .4 l +a 0 .4° 
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where [a)=[A | C | T | G] take the values [0 | 1 1 2 | 3 ] respectively. 

An element of the dot spectrogram is referred to herein as a hixel : h(ij). 

5 The complete spatial randomness (CSR) hypothesis for the dot spectrogram spatial point 

pattern is violated by the microarray oxel layout to maximize posteriori detectability. This means 
that (i) Correct Binding (CB>intensity i.e., number of spatial events (hixels with correct bindings) 
do not follow a Poisson Distribution; and (ii) sampled bindings do not represent independent 
random sample from a uniform distribution on POBC. 

1 0 Compliance to CSR hypothesis would imply constant CB-intensity over the DST, as well 

as no interactions (inhibitions or reinforcement) among the non-zero intensity hixels. Both imply 
an idealized microarray design with positive SNR performance. Thus proclivity for aggregation 
during hybridization is assumed. Computationally, this method assumes that a 1000 X 1000 
microarray (i.e., N = 1,000,000 oxels) is reduced to <O(100) nodes with the preconditioning 

15 method presented below. DST algorithm on preconditioned hixel array, such as variants of 
Dirichlet tessellation would further reduce it to -0(1 0) with no coding. Consequently, this method 
must be applied to ~0( 1 0) ensembles. For ultra-sensitive detection, it is assumed that the ensemble 
basis system is degenerate, i.e., it has (i) a bounded intensity function; (ii) a bounded radial 
distribution function, (iii) and is anisotropic. The method is however valid if hixel intensity 

20 distribution can be approximated by lattice-based processes, Cox processes, Markov point 
processes as well as homogeneous and non-homogeneous Poisson processes. Accommodating the 
complete spatial randomness (CSR) conditions, within a tessellation, to analyze "degenerate, 
contagious distributions" is actually paradoxical to the lineage from QSR and violate temporal 
linearity and contiguity assumptions for amplification of certain spatial effects. 

25 Quantum Expressor Functions 

As noted, the Quantum Expressor Functions (QEF's) generated at step 202 are based upon 
the DNA chip used to generate the dot spectrogram and based upon the mutation set of interest. 
More specifically, as shown in FIG. 4, the QEF is generated based upon the spin Boson basis 
system by first calculating the Hamiltonian for the system at step 402 then, at step 404, calculating 

30 harmonic amplitudes JPJ for the Hamiltonian. Next, at step 406, an order function (OF) is 
generated, then at step 408, entrainment states are measured of the OF of the ground truth. Finally, 
at step 410, the OF of ground truth is modulated to yield the QEF. 

The use of the spin- boson Hamiltonian permits the exploitation of quantum stochastic 
resonance (QSR) phenomena. Classical QSR is an archetypal example of a phenomenon where 

12 



WO 00/52625 PCT/US00/04076 
quantum noise is exploited to drive order in a quantum-mechanical system, as opposed to Gaussian 
noise in a classical system. The appearance of a resonance requires an asymmetry in the energies 
of the two states. A rate equation can be constructed for the system, such that the dynamics can be 
characterized in terms of transition rates <I> + and O- between the two asymmetric quantum 
5 superposition states, and when the drive frequency and the interwell transition rates are much 
slower than the intrawell relaxation rates. The signal to noise ration (SNR) of such a superposition 
system is given by: 

SNR = * W— IT* 

4 1 + exp[f 0 //r B r 0 ][ [jc B rJJ 

10 Where k b is the Boltzmann's constant, and T is the temperature. The sinusoidal ly 

S = € 0 + S e X COSfiJ, / 

modulated asymmetry energy e is given by: 

C(r) = (« + (/ + r|^,/)/t + (/|^,-oo) ) 
In the above expression the power spectrum S(cd) represents the Fourier transform of C(x), 
containing a roughly Lorentzian broadband noise background and 8-function peaks at <a = 0, the 
15 driving frequency co s , and its harmonics. The measured correlation function for the quantum noise 
is given by C(t) = <n +i (t)n +I (t+ x)>, where each n +l = 0 or 1, and the <, indicates an average overt 
over many data points i taken at equal intervals. 

The probability of being in the + quantum state at t after being in the state q„ at is given by 

20 

QSR occurs for an asymmetric well, but not for a symmetric energy well. 

25 A custom Hamiltonian, for use with step 402, which couples the above system to an 

ensemble of harmonic oscillators is given by 
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z z n n 

The dissipative tunneling of the quantum system can be described by the above 
Hamiltonian. In the above expression, e denotes the asymmetric energy, A is the tunneling matrix 
5 element, and Oj are the Pauli spin matrices, and is a harmonic oscillator creation operator with 
frequencies (o n . 

The information about the effects of the environment is contained in the spectral density 

10 However bulk of QSR results to date are limited to either experimental observations in 

esoteric quantum-mechanical systems (e.g., coherent motion of a wave packet in double quantum- 
well (DQW) semiconductor heterostructures) under specific conditions or explanations of 
experiments under carefully controlled conditions. It does not appear that any algorithmic 
implementations or methods have been designed or implemented to enforce nonlinear resonance in 

1 5 spatial, temporal or spatio-temporal phenomena in discrete systems by coupling such systems with 
another mathematical system. 

Thus, as noted above, an important aspect of the exemplary method of the invention is to 
couple the transformed and preconditioned discrete microarray output to a mathematical model for 
a quantum-mechanical dynamical system with specific properties. Specific exemplary parameters 

20 for use in calculating the Hamiltonian are those proposed by A.J. Legett et ah, Reviews of Modern 
Physics, 59, 1, 1987 and A.O. Caldiera and A.LLegett Annals of Physics, 149, 374, 1983. The 
parameters are important only for an offline simulation of this spin Boson system on a digital 
computer. The empirical observables are then collected and used to estimate and compute spectral 
properties, which are actually used by the method. 

25 It is believed that any experimentally or analytical valid parameters for the above system 

will work with the technique because the robustness of the method depends only on the bulk and 
qualitative property shown by this system and the properties of its power spectra. This is an 
important point since the method is actually substantially immune to nuances and specifies of the 
actual driving mathematical system. 

30 The harmonic amplitudes determine the weights of 8 spikes of an averaged spectral power 

density in an asymptotic state S^co). e refers to the coupling strength and P 0 defines the 
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equilibrium state in the absence of driving force. 

The power amplitudes rj m in the mth frequency component of asymptotic state space are 
calculated at step 404 using 
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f7 w (*V) = 4*|/» M (Q,*) 
and the phase shift is given by 

_ \\m{P m (Cl,e )) 



<p m (Q, s ) = tan" , 
* m [Re 

The analytic for the external force is given by 



Pjfl,e) = —L—^h (-imQ, r ) 
y - imQ 7t 



15 The parameters y, e 0 are predetermined and are design specific. Typically, values of 

0.001 and 0.0001 are used for y and e 0 respectively. In the above expression |P TO | 
determine the weights of the 8 spikes of the averaged spectral power density. 

For particular applications, the QEF is designed by matching the power spectral density 
20 (PSD) amplitude of coded mutation signature to that of the spin-boson system described above so 
that stochastic and deterministic time scales match and so that the time scales couple back to noise 
statistics and degree of asymmetry. The method employs a fully automated iterative conjugate 
gradient relaxation method for Spectral matching between asymmetric base system and coded 
mutation signature. The actual determination of the QEF depends on the specifics of 
25 bioelectronics substrates used for actual analysis. The method is however generalizable to all or 
almost all arrayed embodiments. In addition, the method is highly scalable to array dimensions (as 
the offline design trade-space time does not matter to computational complexity). Since the system 
is an overdetermined coupled system, convergence criteria and stability of relaxation method 
directly relates to downstream resonance effectiveness. 

30 
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Order Function 

At step 306, the order function (OF) of ground truth is generated as follows. The order 
function (OF) is for ground truth wherein ground truth represents a state wherein a positive signal 
to noise ratio (SNR) is expected for hixel intensities of selected oxels. 

A general form for the classical Order Function (OF) is 



-ijoke 
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_ 1 

where i = V-f and h k « \d€k{6^e ltdk9 and 0 < G < 1 . 

0 

The following equation 
i 

h k s ^dm{0)e~ lnike is referred to as the Diado integral. 



15 The OF is calculated using an order match which implies the variance of the density 

function for a specific exponential family approaches 0 (within 0.000001 - 0.0001). The free 
energy for the density function is given by 

20 

where % a and £ h represent state vectors. X b Q represents the random observable in 
symmetric bilinear form, and X denotes the characteristic function. 

The detailed form and the specific results for solving the system are given in D.C.Brody 
and L.P. Hughston, Geometry of Quantum Statistical Inference", Physics Review Letters, 77(14), 
25 pp. 285 1 -2855, 1 996. An important aspect of the exemplary method of the invention in using the 
result to derive the QEF which gives unique advantages in the enhancement process. 

Entrainment States 

Thus the OF derivation is based in Diado's theory of multibranch entrainment of coupled 
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nonlinear oscillators, wherein a number of different entrained states co-exist At ground truth 
POBC locations, an a priori measurement of entrainment is performed (step 408). This is done by 
approximate Daido IntegralO with Z* = PSD maxima at regions where DST boundary > desired 
detection threshold. (For more information regarding the Daido Integral, see Physics Review 
5 Letters, 77(7), 1406-1411.) Premultiplier constants are used to ensure that meets the following 
tests: 

- maxima power spectrum density (PSD) matches; and 

- L 2 norm on even harmonics < e where e = NS-MRF Barrier assuming CSR assumptions. 

10 

The notion of entrainment states is exploited, in part, because the method treats the hixel 
dot spectrogram as a special case of coupled nonlinear oscillators in equilibrium. However, due to 
device imperfections, hybridization degradation and other limitations number of entrained binding 
states (i.e., incorrectly hybridized) coexist. A PSD match is desired where the OF is Z-peaked 
1 5 (but single peaked around each MRC-hixel). Absence of a single peak implies perfect or lack of 
hybridization. It also defines the resonance loci for this method (i.e., where maximum SNR 
enhancement) is obtained. 

Ground Truth Modulation 

20 The OF of ground truth is modulated at step 410 to yield the QEF as follows. Under 

controlled calibration, as stated above, maximal SNR enhancement (optimal resonance) is 
achieved when OF yields a single peak. It is a important design point for matching PSD of 
coupling spin Boson system to the synthetic QEF. The specific form of the QEF to be used is the 
generic OF shown above. So the exemplary method exploits two connotations of OF: (a) 

25 parametric form for the QEF (that is closer to the classical form) and (b) as exponential attractor 
for a dissipative system. The two OF's are then recoupled and convolved with the preconditioned 
dot spectrogram (see below). 

The resulting QEF generated via the steps of FIG. 4 is given by: 

30 

QEF MRCi = Hm = - £ ±h a Z,e— 
Typically, only first three PSD peaks are considered for spectral matching. 
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Preferably the QEF is represented digitally using a matrix or array having the same number 
of elements as the dot spectrogram to be analyzed. 

5 Phase Space Resonance Stimulus Using Palm Operators 

Referring again to Fig. 2, if the dot spectrogram generated at step 206 is not a phase-space 
representation of the output of the hybridization chip, then it is desirable to convert the 
preconditioned dot spectrogram (generated at step 208) and the QEF into phase space to facilitate a 
phase space resonance interaction. In other embodiments, though, an amplitude-based resonance 

10 interaction is performed and hence it is not necessary to convert to phase space. In still other 
implementations, other types of resonance interactions may be employed. In the following, it is 
assumed that the dot spectrogram has phase space components and hence the following conversion 
steps are applied to the QEF at step 204 and similar steps are applied to the dot spectrogram 
following image preconditioning. 

15 A phase embedding operator, I\ is applied to the hixels corresponding to the coded base 

mutation set such that hixel values now correspond to angles and not to intensities. These values 
are cyclic, with absolute magnitude of the phase image having no meaning. The relative 
magnitudes are more significant. So if X is a phase image and a is any constant, then X+a (mod 
coIor_map_scale) is a valid descriptor for the phase image. More importantly, the difference 

20 between a max_hixel_intensity and min_hixel_intensity is 1 and not color_map_scale. The 
difference between phase values at two hixels i and j is X(i)-X(j) (mod color_map__scale). 

The phase embedding operator is designed such that transitive closure between any two 
hixels is maintained, i.e., there is an accumulated phase function \\? for which 

25 X(i) « \j/(j) (mod color_map_scale) 

The phase values of \|/ accumulate rather than cycling back to 0. Note that \\t can be 
approximated at least locally by a linear function / with deviation error function err(p) having 
values relatively close to zero for which 

30 

V(p) = /(p) + err(p) 

This function / is the analog of an amplitude image having a constant intensity and err(.) is 
the analog of the deviation from this constant amplitude. It is assumed that err(p) never crosses a 
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phase discontinuity and can be treated like a real-valued function. If values are computed for this 
deviation function err(.), then |err(p)| would specify an intensity image whose local maxima might 
indicate detection of expressed oligonucleotides of interest. 

Each phase angle corresponds to a point on the unit circle in the coordinate lane. This is a 
5 one-to-one mapping in both direction. Let V be the mapping that takes a phase angle into the 
corresponding vector on the unit circle. Let I be the inverse mapping that takes a vector on the unit 
circle back into the corresponding phase angle. V(f(p)) is thus a vector valued function which has 
no associated phase discontinuity. V can be averaged over a region without the problems 
associated with directly averaging phase angles. 

1 0 Given a rectangular neighborhood N(p,m,n) of the hixel p, let v be the vector which is the 

average of all the values of V(f(p)) in this region and let |v| be the magnitude of this vector. 
Clearly |v| must lie between 0 and 1 , inclusive. Moreover, |v| is a measure of the dispersion of the 
phase angles. If |v|=l , then all the vectors in the region N(p,m,n) are equal to v. If |v|=0, then the 
vectors are distributed more or less uniformly about the unit circle and there is no "average" phase 

1 5 angle. If |v|0, then v/|v| is a vector on the unit circle and I(v/|v|) can be defined to be to the 
"average" phase angle ave(f,p,m,n). Clearly, the closer |v| is to 1 , the tighter the clustering of the 
vectors about v/|v| and the greater the certainty that this is a meaningful average value. The 
important characteristic is that this average will tolerate phase with data noise. 

If the line is taken which is perpendicular to the vector v and passes through the point v, 

20 then this line will specify a chord of the unit circle. The arc of the unit circle which corresponds to 
this chord is the region from which the "average" vector comes. The length of this arc is 
2*arccos(|v|). Thus, to ensure that the "average" vector lies in a single quadrant or a 90 degree arc; 

|v| = cos(45 degrees) = 0.7071 

25 

|v| is referred to herein as the average magnitude and the minimum acceptable value for this 
average magnitude is referred to as minimum average magnitude. 

Averaging Phase Input: Let f be a phase image, F the corresponding accumulated phase 
30 function, and L+e=F be as described above. Assume that m and n are odd integers so that 
N(p,m,n) is a symmetric mXn rectangle of hixels with p in the exact center. Since L is linear and 
N(p,m,n) is symmetric about p, the average of L over N(p,m,n) will be L(p). Since e = F-L, the 
average of e over should be N(p,m,n) to be close to 0. Hence, approximately, 
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ave(F,p,m,n) = ave(L+e,p,m,n) 

= ave(L,p,m,n) + ave(e,p,m,n) 
= L(p) + ave(e,p,m,n) 
= Up) 

Since f(p) = F(p) (mod), the average values of these two functions overN(p,m,n) should 
also satisfy this relationship and hence, approximately, 



ave(f,p,m,n) = ave(F,p,m,n) (mod) 
10 =L(p)(mod) 

In other words, it is reasonably expected that ave(f,p,m,n) recovers an approximation to the 
underlying linearity. If this is valid, then 

15 | ave(f,p,m,n) - f(p) | 

is a measure of the deviation of f at p from this underlying linearity and these values form 
an intensity image whose large values and local maxima reasonably correspond to mutations of 
interest. 

20 A possible disadvantage of this approach is that it may be valid only if the phase values 

being averaged are all within a single phase cycle and, preferably, within a fraction of a phase 
cycle (e.g., 120 degrees). If the available phase image is too grainy, it may not be possible to 
average over a large enough rectangle. In this case, it is probably necessary that there be some 
larger scale linear regularity which extends over at least several consecutive phase cycles. The 

25 remaining sections presents alternatives based on this assumption. 



Define the first vertical difference of a phase image f to be 

dv(f,p) = f][i+lj)-f(ij)(mod) 
where p=(ij). Similarly, define the first horizontal difference of a phase image f to be 

dh(f,p) = f(ij+l)-f(ij) (mod) 
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These both define new phase images derived from f and are discrete analogs to the first 
derivative. The range of values for dv and df are also from -1 80 to +1 80. 

Let F be the accumulated phase function corresponding to f and let L+e=F be as described 
above. For the sake of convenience, assume that p is translated to the origin so that p=(0,0). The 
5 linear function L is expressed as 

L(ij) = V*i+H*j + K 

for some constants V, H, and K. If L is a reasonable approximation to F, then it is expected 
1 0 that the first vertical difference of f to be an approximation to the slope V and the first horizontal 
difference of f to be an approximation to the slope H. Hence it is possible to estimate V and H by 
averaging dv and dh over some appropriate rectangle centered at p. Let 

V 1 = ave( dv(f,p), p,m,n) 
15 IT = ave( dh(f,p), p,m,n) 

be these estimates. Note that these are phase averages. Probably m and n should be chosen 
so that m/n is roughly equal to HW\ This ratio would imply that the total phase accumulation 
across the rectangle N(p,m,n) is about the same in both the vertical and horizontal directions. For 
20 example, the rate of change of phase along the i-axis (vertical) may be twice the rate of change 
along the j-axis (horizontal). In this case, a 5x1 0 rectangle of hixels may turn out to be the most 
appropriate neighborhood. 

If a good estimate K' is computed for K, then an approximation to L is: 

25 

LXij^V^i + H^j + K' 

If L* is a good enough approximation to L and F is linear enough, then the accumulated 
phase function F is reconstructed as follows: 
30 e'(p) = f(p)-L'(p)(mod) 

and 



F(p) = L'(p) + e'(p) 
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Note that this may require that the error function e 1 never deviate far enough from 0 to cross 
a phase discontinuity. 

With an accumulated phase function F now explicitly computed on some rectangle 
N(p,m,n), standard numerical analysis techniques are applied, such as linear regression, to find a 
5 function L" which best approximates F over this rectangle. The value 



dev(f,p^Ti,n) = |F(p).L"(p)| 



can be taken as a measure of the deviation at p of the image f from its "expected" value at 
10 p. In particular, this defines a new image whose values are intensities rather than phase angles and 
whose local maxima can be analyzed. 

To compute a good estimate K' for K in order to determine the function L', choose a good 
estimate of K' to be a phase value between -1 80 and +1 80 which minimizes the average value over 
the rectangle N(p,m,n) of either 

15 

|f(p)-L'(p)(mod)| 



or 



20 (f(p)-L'(p)(mod)) 2 

If the accumulated phase function F corresponding to the phase image f is approximately a 
linear function and, as a result, the first differences are approximately constant, then the first 
differences of the first differences should be approximately 0. 
25 Define the second vertical difference of a phase image f to be 

d2v(f,p) - [f(i+l j) - f(i j) (mod)] - [f(ij) - j) (mod)] (mod) 

and the second horizontal difference of a phase image f to be 

30 

d2h(f,p) = [f(i j+1) - f(ij) (mod)] - [f(i j) - f(i j-1) (mod)] (mod) 



This corresponds to an indexing of a one dimensional difference table that looks like the 
following: 
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a0 al a2 a3 a4 a5 a6 
bO bl b2 b3 b4 b5 
cl c2 c3 c4 c5 

5 

where bO=al-aO, bl=a2-al,... and cl=bl-bO, c2=b2-bl,... 

Let F=L+e be define as earlier. If e were identically 0, then F would be linear and both 
second differences d2v and d2h would also be identically 0. If e is not identically 0, the second 

1 0 differences will necessarily also deviate from 0, however, not necessarily in the same place as e. If 
|e| attains a local maximum at a hixel p and this is associated with a very localized deviation from 
0, as would be expected of a signal from a mutation of interest, then d2v and d2h should also 
deviate from 0 fairly nearby to p. 

The deviation in d2v and d2h associated with a deviation in |e| may be spread out over a 

1 5 number of hixels. This will not only make it harder to detect a significant change in d2v and d2h 
but also make it harder to locate the center of this change using a local maxima finding filter. It 
may be possible to compensate for this by averaging d2v and d2h or |d2v| and |d2h| over 
rectangular regions N(p,m,n). This should tend to sharpen the peaks in d2v and d2h. Correlating 
several different images generated using d2v and d2h may also help with center finding since the 

20 correlation algorithm will automatically average nearby local maxima. 

Thus a method for achieving induced coupling for phase space has been described. 
Repeatability of the overall method hinges significantly on the aforementioned coupling as it 
induces analytic redundancy via complex information fusion at constant feed forward 
computational cost. Hence, the coupling is advantageous for at least two reasons: 1) it applies to 

25 any microarray device that may provide a phase or cyclical (modulo) input as opposed to 
amplitude input; and 2) it is important for use with active hybridization devices which will have an 
element of built in control that will have a phase representation. The coupling is actually 
introduced into the exemplary method in three possible ways. 

(a) if the entire resonance computation proceeds in the phase space as opposed to amplitude 
30 space. Then it is necessary to transform QEF and post-hybridization microarray output to phase 

space. The coupling material provides the transformations to do so. 

(b) if the microarray output is already in the phase space then it will only be applied to 
convert the QEF. 

(c) if different oxels are designed to solve for different aspects of the analysis, the coupling 
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permits a combination of outputs to reach a final conclusion. 

Entrained states are employed, in part, to precisely compensate for that point. So already in 
the method synthetic decomposition and coupling are provided. But it is used to reject spurious 
candidates as opposed to "light up" more oxels with potential match. 
5 Again the step of converting to phase space may be optional depending upon the 

implementation and is applied, if at all, after preconditioning and/or before introducing the QEF to 
the resonant interaction step. 

Also, a major limiting restriction in QSR that is avoided by the exemplary method pertains 
to matching the stochastic and deterministic time scales in "domain system" and the external 
1 0 coupling asymmetric dynamical system, since this has limited applicability to continuous data. By 
replacing the time scale match requirement with ensemble spatial statistics (structure of 
oligonucleotide itself) expressed via the aforementioned generators an entirely new class of 
analysis for discrete systems is enabled By deriving the condition that labeled oligonucleotide 
chains need to satisfy to maximize QEF generator capacity, a properly constructed and conditioned 
15 QEF will be able to extract and enhance an entire class of mutations over and beyond single 
mutation detectability. 

The resulting phase-space representation of the resonance stimulus generated at step 204 is 
given by: 

20 ^n^c stimulus = (D2h(QEF MRC { ) <8> mod „ D2v(QEF MRC j) for all QEF subarray elements. 

As with the QEF, the resonance stimulus pattern is preferably represented digitally using a 
matrix or array having the same number of elements as the dot spectrogram to be analyzed. 

25 Generation of The Dot Spectrogram 

Continuing with the detailed description of the steps of FIG. 2, as noted a dot spectrogram 

is generated at step 206 for a sample from an N by M DNA chip array wherein an element of the 

array is an "oxel": o(ij). A 6-cr manufacturing process accuracy in microarray design is assumed. 

Each array cell amplitude is given by *(i j) for i: 1 to N, and j: 1 to M. Let \j/(i j) denote the a 
30 priori known oligonucleotide given by [ct^x,...ct k ], where a = {A,C,T,G} base associated with 

each array eel I [a,b] where 1 0 < k < 25 . The complimentary strand, derived from unknown sample 

is denoted by ¥(/,/) . 
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The post-hybridization microarray is treated mathematically using the machinery of 
equations with aftereffect. Each hixel given by 0(i j) is represented as a cluster of dynamical 
systems of potentially [CB] correctly bound, [UB] unbound, [PB] partially bound and [IB] 
incorrectly bound. Thus [CB] ^ + [UB] ^ + [PB] <d (U) + [IB] ^ « To^within 0.0001%. 
5 The model analytic for the estimated fluorescence activity T in a hixel h^ is given by 




t tp (hybridization start time). 



The fluorescence stabilization section is given by 

10 

And the rate of fluorescence stabilization is expressed via 
15 4* = 4* (t+6) 9 S<0; 

While , A(f >} = ^ ^ (/ + S) , £ < 0; the system is assumed 
to be memoryless, i.e., 0 is ergodic. 

20 The resulting dot spectrogram generated at step 206 is given by: 



25 

Preconditioning of The Dot Spectrogram 

With reference to FIG. 5, the dot spectrogram <t>(i j) is preconditioned by performing the 
following steps. First, at step 502, the dot spectrogram is refocused to yield a refocused dot 
30 spectrogram. Then, at step 504, a cross-correlation convolution operation is performed to yield a 
correlated refocused dot spectrogram. A local maxima filter S is then applied at step 506 to the 
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correlated refocused dot spectrogram to yield a maximized dot spectrogram. The maximized dot 
spectrogram is re-scaled at step 508 to yield a uniformly re-scaled dot spectrogram. The uniformly 
re-scaled dot spectrogram is then further re-scaled at step 510 by amplifying local edge hixel 
boundaries of the uniformly re-scaled dot spectrogram to yield a globally re-scaled dot 
5 spectrogram denoted by &(JJ). Also, amplitude wanderings within the globally re-scaled dot 
spectrogram are estimated at step 512 for use downstream. 

The steps of FIG. 5 result in a dot spectrogram tessellation (DST) operator. Purpose of 
DST Operator is to determine idealized ensemble boundaries for forcing downstream resonant 
action. Typically, in signal processing applications, high pass or band pass spatial filtering is 

1 0 implemented to enhance SNR in dot spectrogram matrix. Alternate methods apply a combination 
of Laplacian or other edge detection filters apply to enhance signal from arrays cells with a higher 
hybridization concentration from those of the adjacent cells. These SNR enhancement methods 
however work only with positive or zero-SNR. Since SNR in general is negative in our case (ultra- 
low target DNA concentrations), these methods in effect amplify noise or further blur the hixel 

1 5 boundaries. 

With the exemplary method, a new technique is used to amplify highly local morphological 
variations in the overall dot spectrogram matrix. A filter £ is applied to the normalized amplitude 
dot spectrogram matrix <I>, for all combinations ±i±j±k where k is typically ranges from 0 to 
2. 

20 

Refocusing The Dot Spectrogram 

Refocusing of the dot spectrogram at step 502 to yield a refocused dot spectrogram is 
performed by determining a locally averaged amplitude sub-array represented by Q, then for 
each value of i, j and k, where k ranges from 0 to 2: determining a local standard deviation a 
25 over a (2k+l)X(2k+l) hixel neighborhood centered at (ij) and applying a filter ^ to the dot 
spectrogram O, for all combinations ±d±j±k. 

a is calculated by 



30 ctZ =Jz z(*>(ij)-n z y 

\i-2 J-i 
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The calculation of CI is independent of the oxel layout. From a preconditioning 
effectiveness standpoint, a best case design corresponds to a completely random oxel layout in 

terms of O -value for adjacent oxels. The worst case corresponds to O -value separation of 1 
among adjacent oxels. The oxel, o(ij) 2 , centered at (ij) comprises of complementary 
5 oligonucleotides, corresponding to a mutation of interest over the set Z. Loci, r, for averaging 
amplitude ranges from ±5 oxels to ±k oxels depending on ( 0-value mod 4*) separation. Qf then 

is computed using Qf = £ zl 'j? 2 v z - Wh en the averaging loci for multiple mutations of 

interest overlaps, then composite loci, r' is bounded by a rectangle whose top left hand and bottom 
right hand corner coordinates are given by [u-k-1, j^+k+1] and [i mix +k+l, j min -k-l] where i min , 
10 'max» jmin anc * jmw correspond to the ordinate and abscissa for the oxel mutations with overlapping 
loci. 



The preceding expression captures local ensemble deviation from dot spectrogram average. 
The design parameter k is computed offline based on specific bio-molecular signatures of interest. 
The suffix / denotes iterative index for a cross-correlation convolution operation to be applied after 
20 global refocusing. 

A uniform dot spectrogram rescaling can be achieved by applying a . 3(1 j) where a can 
be either a constant or a functional. This operator selectively enhances those hixels and ensemble 
boundaries whose intensity exceeds the local average by more than k prespecified standard 
deviations. 

25 k is a predetermined constant computed offline based on microarray fluorescence or 

chemiluminiscence sensitivity. 



£ given by: 



15 



<l>(iJ)i = ^(ij)-ic-[n(0(ij),k)/a(0(ij),k)] V ij 




30 



The suffix / denotes an iterative index for a cross-correlation convolution operation to be 
applied after global refocusing. 
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Cross-correlation of the Dot Spectrogram 

The cross-correlation convolution operation of step 504 to yield a correlated refocused dot 
spectrogram is performed by convolving the refocused dot spectrogram with an apriori-chosen 
restricted field obtained using a strongly dissipative dynamical system. Dissipative dynamical 
5 systems are those which define a forward regularizing flow in an adequate phase space containing 
an absorbing set. An absorbing set is a bounded set that attracts all bounded solutions in finite time 
at an exponential rate. Since we exploit a strongly dissipative system, the absorbing set is required 
to be unique compact set that is both positively and negatively invariant under the flow, it attracts 
all flows. 

10 Conceptually, the post-hybridization hixel denotes a projection of the flow field that is 

absorbed by the coupling system. Since we couple with an absorbing set, this stage yields 
significant SNR enhancement. 

In general, the convolving field can be constructed using Kuramoto-Sivashinsky equation, 
2D Navier-Stokes equation or some forms of Reaction-Diffusion equations. The dot spectrogram 

1 5 subarray around the oxel detecting mutation of interest can be cross-correlated with any dissipative 
dynamical system. In summary this step exploits proven classes of mathematical system or 
exponential attractors. By definition exponential attractor is an exponentially contracting compact 
set with finite fractal dimension that is invariant under the forward flow. 

In the exemplary implementation, the following is performed to achieve the cross- 

20 correlation. A synthetic system of the form is assumed, whose evolution equation is given by 

u t +Au + R(u) = 0 
u(0) = u o 

25 where A is a positive self-adjoint operator with compact inverse and nonlinear term R 

satisfies the order match condition in step 306 in Figure 3. 

The functional formulation of 2D Navier-Stokes equations for incompressible fluid flow is 
chosen as the coupling system 

30 
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It is given by 

m, + vAu + B(u 9 u) = /, 
u(0) = u o 

where A = -P H A is the Stokes operator, B(u,u) stands for the non-linear term (u . Au) 
projected to the underlying Hilbert space H, f is the volume force projected to the same Hilbert 
space and v is the viscosity term. To meet the incompressibility condition the following Hilbert 
space is assumed to set the initial value problem 

H = | u e ]} (e) 2 : div u = 0 , Ju(x)dx = 0 , u,| ^ = u,| „ =0 , J = 1,2 
and 



10 



15 



where Q is the square [0,L] X [0,L]. The mathematical domain of the Stokes operator is 
given by 

D(A) = H 2 (Q) 2 nV 

Gaterkin Approximation can be used to approximate the exponential attractor for the 

system 



The above system with periodic boundary conditions admits an exponential fractal attractor 
20 M in B whose dimensions can be estimated using 

dF{M) < cG 2 (logG 4 vX, +1) 

where G is the Grashoff number given by !j~ and c is a constant that depends on the 

V Ay 

shape factor. X t denotes the positive eigenvalue of A. 
25 The rate of convergence of this system can be computed as well. 

The exponential attractor is then coupled with the post-hybridization dot spectrogram 
subarray 
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The exponential attractor is discretized over the grid that corresponds to the refocused 
amplitude subarray associated with a mutation and estimated above. 
The actual convolution of the two systems is then given by 

5 

-rZS«P(- s Rj)*«('.-0) 
r * j 

This process is computed for each oxel associated with the loci of a specific mutation of 
interest. 

10 

Maximizing the Dot Spectrogram 

Step 506 for applying a local maxima filter 3 at to the correlated refocused dot spectrogram 
to yield a maximized dot spectrogram is then implemented as follows. The local maxima filter is 
defined by 

15 

wherein 3(1 j) = 1 if 4>(i j)' is a local maxima (or defines a local maxima ensemble with 
several neighboring hixels with same amplitude intensity); and wherein S(I j) = 0 otherwise. 

20 

Re-scaling the Dot Spectrogram 

The maximized dot spectrogram is rescaled at step 508 to yield a uniformly re-scaled dot 
spectrogram by applying an operator a . 8(1 j) where a is either a predetermined constant or a 
predetermined functional. 

25 

a is set so as to uniformly scale the convolved dot spectrogram to l/( l+exp^MINn^)) to 
l/(l-texp(MAX INT /2)). Where [MIN 1Nrr , MAX 1NT ] refers to the dynamic range of the microarray. 

The uniformly re-scaled dot spectrogram is itself then rescaled at step 410 by amplifying 
30 local edge hixel boundaries of the uniformly re-scaled dot spectrogram to yield a globally re-scaled 
dot spectrogram. This is achieved by 1 ) determining the zero mean amplitude for the uniformly re- 
scaled dot spectrogram; 2) applying a logarithmic rescaling function p around the zero mean 
amplitude; and 3) merging the local maxima into a single local maximum halfway in between. 
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The logarithmic rescaling function is generated by generating an expansion sequence of 
nonnegative numbers and by generating an expanded dot spectrogram tessellation for <X>. 



The expansion sequence is generated as follows: 
(P}= Po» &\> &2> ». 

which is strictly decreasing until it reaches zero and thereafter is equal to zero by using: 



The expanded dot spectrogram tessellation for 3> (which is represented by DST(O)) is 



wherein D(<k, <J> k ) is a shortest possible discretized fluorescence amplitude separation 
between a pair of hixels <J>j, <|> k wherein <J> k is a local maximum. 

The local maxima are merged into a single local maximum half way in between for 
downstream hixel-to-ensemble and ensemble-to-ensemble operations on hixel clusters using the 
sequence 

max{ X(c)- K(m),0 }; m = 0,1,2,... 
where c is some positive integer constant. 

Additional details are as follows, for any pair of dot spectrogram hixels, <J>j and fy, let the 
discretized fluorescence ampl itude separation between them be denoted by D^, <j>j). So, D(4>j, <j>j) = 
0 if and only if 



I4>s • 4>jl ^ C where C = Max( p(0>(i j))) - Min( p(<D(i j)))/N and D(<|>„ <t>j)=l if and only if ^ 




for 7i > 0 
for/?<0 



generated using: 



DSTCO) -*(£>(*,, + k » 
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and fa are adjacent to each other on the normalized scale, but not identical, i.e., \fa . fa\ ^ 2C. 

Let 4>k be a local maxima image for some other dot spectrogram realization o(i 9 jfi * Define D(4> b 
<J> k ) to be the shortest distance possible from hixel fa to some hixel (|> k which is a local maximum. 
Clearly, DOh, = 0 if and only if fa, is a local maximum and D(fa> fa,) = 1 if and only if <J) k is not a 
5 local maximum but is adjacent to some local maximum. 



The an expansion sequence is thus defined to be any sequence of nonnegative numbers { gy } 
= Po> &>if &2> — (as mentioned above) which is strictly decreasing until it reaches zero and 
thereafter is equal to zero. For example, { p } = 1 5, 1 1, 5, 2, 1 , 0, 0, 0, 0,... If O is a local maxima 
1 0 DST and { p } is an expansion sequence, then the expanded DST(4>) is defined as 



DST(<D)=p(D((l) i , * k )) 



K: Z -»Z on integers is useful for defining expansion sequences. A key property of this 
1 5 sequence is that given only the value of the sum s=p+q, where p and q are any numbers from this 
sequence, then we can always determine the values of the two original summands p and q. For 
example, choose the expansion sequence to be: 

31, 15, 7, 3, 1, 0, 0, 0, 0,... (generated by X(5-i): for i = 0,1,2, ...) 

20 

and let DST A and DST B be two DST realizations "expanded" using this sequence. This 
unique summation property ensures that a threshold value can be chosen for the DST summation 
DST A +DST B so that the resulting realization retains only those local maxima in the two images 
which are within bounded intensity of n of each other. This expansion sequence preserves the 
25 individual local maxima from each of the two realizations. 

For downstream hixel-to-ensemble and ensemble-to-ensemble operations on hixel clusters 
the local maxima are merged into a single local maximum half way in between using the sequence 
using 

max{ N(c)-tt(m), 0}; m = 0,1,2,... 

30 

where c is some positive integer constant. 

In the alternative, a Dirichlet tessellation operator or a Delaunay tri angulation operator are 
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then applied to perform gradient refocusing rather than steps 502 - 510. 



5 Estimation of Amplitude Wanderings 

Amplitude wanderings are estimated at step 412 within the enhanced dot spectrogram. The 
estimate is performed by applying a Palm Distribution operator to the globally re-scaled dot 
spectrogram to capture amplitude wanderings and transitions at element, neighboring pair and 
local ensemble levels. The application of the Palm Distribution operator generates bounds that are 
1 0 used to accommodate degradation of hybridization over time. Essentially the estimate exploits the 
use of generator functions to capture stochastic variability in hybridization binding efficacy and 
draws upon results in stochastic integral geometry and geometric probability theory. 

Geometric measures are constructed to estimate and bound the amplitude wanderings to 
facilitate detection. In particular we seek a measure for each mutation-recognizer centered (MRC-) 
1 5 hixel that is invariant to local degradation. Measure which can be expressed by multiple integrals 
of the form 

m(Z)=[f(z)dz 

20 where Z denotes the set of mutations of interest. In other words, we determine the function 

f(z) under the condition that m(z) should be invariant with respect to all dispersions Also, up to 
a constant factor, this measure is the only one which is invariant under a group of motions in a 
plane. In principle, we derive deterministic analytical transformations on each MRC-hixel., that 
map error-elliptic dispersion bound defined on <R 2 (the two dimension Euclidean space - i.e., oxel 

25 layout) onto measures defined on 9t. The dispersion bound is given by 

Log 4 (0 (/J) 

Such a representation uniqueness facilitates the rapid decimation of the search space. It is 
30 implemented by instantiating a filter constructed using measure-theoretic arguments. The 
transformation under consideration has its theoretical basis in the Palm Distribution Theory for 
point processes in Euclidean spaces, as well as in a new treatment in the problem of probabilistic 
description of MRC-hixel dispersion generated by a geometrical processes. Latter is reduced to a 
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calculation of intensities of point processes. Recall that a point process in some product space E X 
F is a collection of random realizations of that space represented as {(e;, Q, | e, e E, fj € F }. 

The Palm distribution, n of a translation (TJ invariant, finite intensity point process in 
is defined to the conditional distribution of the process. Its importance is rooted in the fact that it 
5 provides a complete probabilistic description of a geometrical process. 

In the general form, the Palm distribution can be expressed in terms of a Lebesgue 
factorization of the form 

E/>N*= Ai^xn 

10 

Where n and A completely and uniquely determine the source distribution P of the 
translation invariant point process. Also, E/> N* denotes the first moment measure of the point 
process and is a probability measure. 

Thus we need to determine n and A which can uniquely encode the dispersion and 
15 amplitude wandering associated with the MRC-hixel. This is achieved by solving a set of 
equations involving Palm Distribution for each hybridization (i.e., mutation of interest). Each 
hybridization is treated as a manifestation of a stochastic point process in 9? 2 

In orderto determine II and A we have implemented the following measure-theoretic filter: 

20 Determination of A 

using integral formulae constructed using the marginal density functions for the point 
spread associated with MRC-hixel(i j) 

The oligonucleotide density per oxel p m{fJ) , PCR amplification protocol (a m ), fluorescence 
25 binding efficiency (r\ m ) and imaging performance (xa m ) provide the continuous probability density 
function for amplitude wandering in the m-th MRC-hixel of interest. Let this distribution be 
given by (p m(IJ) , <J m> r lmy ar m ) . 



The method requires a preset binding dispersion limit to be provided to compute A. The 
30 second moment to the function 

at SNR = 0 condition is used to provide the bound. 
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Determination of n 



Obtained by solving the inverse problem 
IT = 0 + P 

where 



where x, and x 2 represent the normalized hybridization dispersion limits. These number are 
empirically plugged in. We choose 0.1 and 0.7 respectively to signify loss of 10% - 70% 
10 hybridization. 

Also , & denotes the distribution of known point process. We use the form 
l/(l+exp( £>(...))) to represent it. 

Ultimately, the final preconditioned (or enhanced) dot spectrogram generated by the step of 
15 FIG. 5 is represented by: 



20 Resonant Interaction 

Referring again to Fig. 2, at step 210, the resonant interaction between the QEF and the 
preconditioned dot spectrogram is performed until a pre-selected degree of convergence is 

ecu = - dV W + + QN(t) 
du 

achieved. The resonance interaction proceeds by iteratively calculating a resonance equation until 
convergence is achieved. 



25 



More specifically, the following resonance equation is iteratively calculated. 



ecu = - dV W + + QN(t) 
du 
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In the above equation V(u) is the actual precondition, refocused MRC-hixel subarray. So a 
stable equilibrium state (microarray) is transformed and modulated with the QEF, i.e., 
5 mathematically destabilize it to achieve a nonlinear resonance point. 

In the forgoing, "a" is a constant « fluorescence decoherence timescale(s). State variable 
u corresponds to nonstationary Markov random field (NS-MRF). Y(t) = Asin (wt + 1) where t is a 
small random phase factor. "A" is a gain function (control parameter that influences convergence 
rates), "t" denotes integration timestep. QN(t) corresponds to the log (PSD maxima) of the ground 
1 0 truth order function summed over all regions of interest. More precisely, QN(t) is reverberation 
projection (at some instant t prior to resonant convergence) for the coupled Ofs discussed above in 
relation to step 204. 

The Palm generators are used as additive correction terms to the potential gradient to 
compensate for uncertainty and post hybridization decay. The actual dynamics is given by V(u, n, 
15 X). 

The resonant iteration is terminated when 



20 
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Log 



£/(/)- t/ f 



£1 



or when iteration counter t exceeds 10 3 (digital approximation to analog dynamics). 

In this method the output readout from the resonance interaction is trivialized to deliver 
point-of-care (POC) diagnostics systems. After exponential SNR enhancement (i.e., resonance) 
only those k MRC-hixels will have an amplitude which are actually present in the unknown target 
sample. 



In the event no resonance is achieved during step 210, a new QEF (corresponding to a new 
mutation of interest) is selected and the method reinitializes to the original DS as computed in step 
508 in FIG. 5. This technique, referred to herein as a "software programmable QEF flush control 
methodology", of reinitializing the dot spectrogram and using a new QEF leads to a new cycle 
30 through a signal sorter for rapid single- / multi-point gene/mutation sorting. It is amenable to RF, 
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electronic as well as optoelectronic QEF loading to microelectronics computing backplane, 
microarray readout and analysis backplane itself optoelectronically or electronically bonded to 
bioelectronic substrate. Furthermore the method can be implemented on offline miniature custom 
VLSI / palmtop / desktop setup. The QEF is also implementable in electronic, optoelectronic, 
5 bionic and ionic COTS/custom device physics. In the event the chip is designed such that 
complementary oligonucleotides for several mutations are spatially well spaced out that multiple 
QEF can be introduced to stimulate the system at the same time, then a logical "OR"ing of the 
resonance output is used to reach the detection decision. 

Also, note that a chip fabricated to implement the exemplary method can work in two 

10 ways: a) seek all mutations of interest simultaneously; or b) seek all mutations of interest serially. 
An advantage of (a) over (b) is computational speedup. An advantage of (b) over (a) is 
serendipity, i.e., in method (a) only those mutations are resonantly amplified that are detection 
candidates. Everything else is likely suppressed or decimated. In (b) multiple resonant outputs can 
be accepted. By accepting multiple peaks for the OF the method can actually accept 2 nd order, 3 rf 

1 5 order entrained states, where order implies hamming distance to the mutation of interest in terms 
of base pair labels and locations. This can be used to theoretically accept entire families of derived 
mutations. 

Preferably, the resonance interaction is performed digitally by applying a matrix 
representative of the resonance equation to a matrix representative of the resonance stimulus in 
20 combination with a matrix representative of the dot spectrogram. 

The final result of step 210 is a set of hixel locations wherein resonance has occurred 
identified by row and column number (i j). 

Resonance Output Interpretation to Identify Diseases 
25 Once resonances peak are observed in specific hixel locations of the preconditioned dot 

spectrogram at step 210, the hixel addresses (k,l) of those locations are mapped at step 212 into the 
oligonucleotide table mentioned above (which contains the oligonucleotide sequences associated 
with hixel locations) to thereby identify the mutations, if any, present in the sample being 
analyzed. This is a simple table look-up resulting in a direct readout of the mutations. For a 
30 custom POC diagnostic sensor only those hixels which relate to mutations or expressed products of 
interest are stored in the table. 

At step 214, the mutations identified at step 212 are then mapped into the mutations table 
mentioned above (which contains the diseases associated with the mutations) to thereby identify 
the diseases, if any, present in the sample being analyzed. This is also a simple table look-up 
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resulting in a direct readout. Again, for a custom POC diagnostic sensor only those mutations 
which relate to diseases of interest are stored in the table. 

More specifically, at steps 212 and 214, the resonant output interaction is interpreted to 
5 yield a set of confirmed mutations as follows: 

MRC-hixels z = o(h(k,l)) where o is some hashing function or table look-up. The step will 
readout the oligonucleotide sequence from a table that has encoded the microarray. Note that this 
mapping also operates to map the confirmed mutations to known diseases associated with the pre- 

1 0 selected set of mutations of interest to identify diseases, if any, indicated by the DNA sample (step 
214). Note also that no probabilistic inferencing, exploitation of learning or nonlinear mapping is 
required to interpret the resonance output. Rather interpretation is rendered very straightforward 
thereby requiring only a low-cost hardware implementation with simple software to implement 
steps 212 and 214. So an important aspect of the exemplary method is to accurately and robustly 

15 detect specific oligonucleotide sequences in the target sample. Subsequent association to 
understood genomics pathways is trivialized. 

Diagnostic Confirmation 

Finally, the diagnosis generated by step 214 is confirmed at step 216 by taking the 
20 identified diseases and solving in reverse for the associated QEF and then comparing the 
associated QEF's with ones expected for the mutations associated with the identified diseases to 
verify correspondence and, if correspondence is not found, then a new set of mutations of interest 
are selected and all steps repeated. 

In other words, this step maps detected mutations and expressed genes to a diagnostic 
25 assessment This is a probabilistic or deterministic step, depending upon the genomics of the 
specific disease. It is represented as 

Diagnosis P FA ^, m ^mutation, mutation^ mutation,) 

Where $ is the specific diagnostic model associated with a medical condition of interest. 
30 Step 216 is a double check mechanism adopted to confirm multi -factorial diseases where 

the biochip encodes complex genomics. 



Alternative Embodiments 

What has been described thus far is an implementation whereby a phase-space resonance 
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interaction is performed. This implementation, which may also be referred to as a "mixed-mode 
phase shifted mode" is particularly effective for automatically extracting an entire class of 
mutations that may be manifested in a hybridized element. In general, the mixed-mode provides 
polymorphism in induced couplings for QEF design which delivers repeatability on analysis 
5 whereby mutation signatures of interest are simultaneously coupled to many base "dynamical 
systems" with a single phase-embedding operator. Other resonance coupling interactions may be 
exploited as well. Other examples of couplings, other than phase-based, are "additive coupling 
mode" which provides further SNR enhancement and a "shunted input multiplicative coupling 
mode" which amplifies noise-to-noise couplings and leads to derivation of better readout threshold 
10 for diagnostics decision making. Also, a combination of different resonance interactions can be 
exploited. 

Details regarding a related implementation may be found in co-pending U.S. Patent 
Application 09/253,792 entitled "Method and Apparatus for Interpreting Hybridized Bioelectronic 
DNA Microarray Patterns Using Self Scaling Convergent Reverberant Dynamics", is incorporated 

1 5 by reference herein. Details regarding an implementation directed to measuring viral loads may be 
found in co-pending U.S. Patent Application 09/253,791 entitled "Method and Apparatus for 
Exponentially Convergent Therapy Effectiveness Monitoring Using DNA Microarray Based Viral 
Load Measurements", is also incorporated by reference herein. 

The exemplary embodiments have been primarily described with reference to flow charts 

20 illustrating pertinent features of the embodiments. Each method step may also represent a hardware 
or software component for performing the corresponding step. It should be appreciated that not all 
components or method steps of a complete implementation of a practical system are necessarily 
illustrated or described in detail. Rather, only those components or method steps necessary for a 
thorough understanding of the invention have been illustrated and described in detail. Actual 

25 implementations may utilize more steps or components or fewer steps or components. 

The description of the exemplary embodiments is provided to enable any person skilled in 
the art to make or use the present invention. Various modifications to these embodiments will be 
readily apparent to those skilled in the art and the generic principles defined herein may be applied 
to other embodiments without the use of the inventive faculty. Thus, the invention is not intended 

30 to be limited to the embodiments shown herein but is to be accorded the widest scope consistent 
with the principles and novel features disclosed herein. 
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CLAIMS 

1 . A method for analyzing an output pattern of a biochip to identify mutations, if any, 
present in a biological sample applied to the biochip, said method comprising the steps of: 

generating a resonance pattern representative of resonances between a stimulus pattern 
5 associated with a set of known mutations and the output pattern; and 

interpreting the resonance pattern to a yield a set of confirmed mutations by comparing 
resonances found therein with predetermined resonances expected for the set of mutations. 

2. The method of claim 1 wherein the output pattern is derived from a biochip having 
1 0 a plurality of cells each having a set of substantially identical immobilized oligonucleotides with 

the oligonucleotides of each cell being substantially unique from the oligonucleotides of every 
other cell. 

3. The method of claim 1 wherein the step of generating a resonance pattern includes 
1 5 the step of iteratively processing the output pattern by performing a convergent reverberation to 

yield a resonance pattern representative of resonances between a predetermined set of selected 
Quantum Expressor Functions and the output pattern until a predetermined degree of convergence 
is achieved between the resonances found in the resonance pattern and resonances expected for the 
set of mutations. 

20 

4. The method of claim 2 further including the steps of: 
selecting a sub-set of the mutations for analysis; 

generating preconditioner transforms based upon the selected sub-set of mutations by 
inverse Dirichlet tessellation on the biochip structure; and 
25 selecting a sub-set of nonlinear Quantum Expressor Functions from a set of predetermined 

nonlinear Quantum Expressor Functions based upon the selected sub-set of mutations. 

5. The method of claim 1 further including the step of modifying the output pattern to 
generate a modified output pattern by performing one or more of the steps of: 

30 refocusing the output pattern to yield a re-focused output pattern; 

cross-correlating the re-focussed output pattern; 

applying a local maxima filter to the correlated re-focussed output pattern to yield a 
maximized output pattern; 

re-scaling the maximized output pattern to yield a uniformly re-scaled output pattern; and 
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re-scaling the uniformly re-scaled output pattern to amplifying local boundaries therein to 
yield a globally re-scaled output pattern. 



6. The method of claim 5 wherein the step of re-scaling the maximized output pattern 
5 is performed by applying an operator to selectively enhance those hixels and ensemble boundaries 

whose intensity exceeds the local average by more than a prespecified number of standard 
deviations. 

7. The method of claim 5 wherein the step of rescaling the uniformly re-scaled output 
10 pattern is performed by applying a logarithmic rescaling function around zero mean amplitudes 

within the output partem to amplify local edge boundaries of the output pattern. 

8. The method of claim 1 further including the step of 

estimating amplitude wanderings within the globally re-scaled output pattern. 

15 

9. The method of claim 8 wherein the step of estimating amplitude wanderings is 
performed by applying spectrogram ensemble invariance property generators constructed using a 
Palm distribution. 

20 10. The method of claim 1 wherein the step of interpreting the resonance pattern to a 

yield a set of confirmed mutations by comparing resonances found therein with predetermined 
resonances expected for the set of mutations includes the step of 

mapping the confirmed mutations to known diseases associated with the pre-selected set of 
known mutations to identify diseases, if any, indicated by the biological sample. 

25 

11. The method of claim 1 further including the step of performing a diagnostic 
confirmation by taking the identified diseases and solving in reverse for the associated Quantum 
Expressor Functions and then comparing those Quantum Expressor Functions with ones expected 
for the mutations associated with the identified disease to verify correspondence and, if no 

30 correspondence is found, selecting a new sub-set of known mutations and repeating all steps. 

12. The method of step 1 further including the initial steps of generating the set of 
nonlinear Quantum Expressor Functions by: 

inputting a set of mutation signatures representative of the pre-selected set of known 
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mutations; 

inputting a representation of a biochip oligonucleotide pattern layout for the biochip from 
which the output pattern is generated; 

generating a set of resonant interaction parameters representative of mutation pattern 
5 interactions between elements of the biochip including interactions from a group including 
element-to-element interactions, element-to-ensemble interactions, ensemble-to-element 
interactions, and ensemble-to-ensemble interactions; and 

generating the set of nonlinear Quantum Expressor Functions from the set 
of resonant interaction patterns. 



13. The method of claim 12 wherein the step of generating the set of nonlinear 
Quantum Expressor Functions comprises the steps of matching a power spectral density (PSD) 
amplitude of a coded mutation signature, corresponding to the pre-selected mutation set of interest, 
to that of Hamiltonian system so that stochastic and deterministic time scales match, and the time 

1 5 scales couple back to noise statistics and degree of asymmetry. 

14. The method of step 1 further including the initial steps of generating the set of 
nonlinear Quantum Expressor Functions by: 

calculating values representative of a pre-selected Hamiltonian function; 
20 calculating harmonic amplitudes for the Hamiltonian function; 



15. The method of step 14 wherein the Hamiltonian is a Spin-Boson Hamiltonian. 

1 6. The method of claim 1 4 wherein the step of generating the order function includes 
the step of approximating a Daido Integral. 



10 



generating an order function from the Hamiltonian; function 
measuring entrainment states of the order function; and 

modulating the order function using the entrainment states to yield the Quantum 



Expressor Function. 



25 



30 



1 7. The method of claim 1 wherein the output pattern is a quantized output pattern. 



1 8. The method of claim 1 7 wherein the quantized output pattern is a dot spectrogram. 
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19. The method of claim 1 wherein the biological sample selected from a group 
consisting of a DNA, RNA, protein, peptide-nucleic acid (PNA) and targeted nucleic amplification 
(TNA) samples. 

5 20. A method for analyzing a biological sample to identify mutations, if any, present in 

the sample from a pre-selected set of known mutations, said method comprising the steps of: 

applying the sample to a biochip to generate an output pattern representative of quantized 
hybridization activity of oligonucleotides in the sample; 
1 0 selecting a sub-set of the mutations for analysis; 

generating preconditioner transforms based upon the selected sub-set of mutations; 
selecting a sub-set of nonlinear Quantum Expressor Functions from a set of predetermined 
nonlinear Quantum Expressor Functions based upon the selected sub-set of mutations; 

preconditioning the output pattern in accordance with the preconditioner transforms to 
1 5 generate a modified output pattern; 

iteratively processing the modified output pattern by performing a convergent reverberation 
to yield a resonance pattern representative of resonances between the selected set of Quantum 
Expressor Functions and the modified output pattern until a pre-selected degree of convergence is 
achieved between resonances of the resonance pattern and mutations of the selected sub-set of 
20 mutations; and 

interpreting the resonance pattern to a yield a set of confirmed mutations by comparing the 
resonances therein with resonances expected for the mutations of the selected set of mutations. 

21. The method of claim 20 further including the step of mapping the confirmed 
25 mutations to known diseases associated with the pre-selected set of known mutations to identify 

diseases, if any, indicated by the biological sample. 

22. A method of generating a set of nonlinear Quantum Expressor Functions 
comprising the steps of: 

30 inputting a set of mutation signatures representative of the pre-selected set of known 

mutations; 

inputting a representation of a biochip; 

generating a set of resonant interaction parameters representative of mutation pattern 
interactions between elements of the biochip including interactions from a group including 
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element-to-element interactions, element-to-ensemble interactions, ensemble-to-element 
interactions, and ensemble-to-ensemble interactions; and . 

generating the set of nonlinear Quantum Expressor Functions from the set of resonant 
interaction patterns. 

5 

23. A system for analyzing an output pattern of a biochip to identify mutations, if any, 
present in a biological sample applied to the biochip, said system comprising: 

a resonance pattern generator operative to generate a resonance pattern representative of 
resonances between a stimulus pattern associated with a set of known mutations and the output 
10 pattern; and 

a resonance pattern interpretation unit operative to interpret the resonance pattern to a yield 
a set of confirmed mutations by comparing resonances found therein with predetermined 
resonances expected for the selected set of mutations. 

15, 24. The system of claim 23 wherein the output pattern is derived from a biochip having 

a plurality of cells each having a set of substantially identical immobilized oligonucleotides with 
the oligonucleotides of each cell being substantially unique from the oligonucleotides of every 
other cell. 

20 25. The system of claim 23 wherein resonance pattern generator includes a resonant 

interaction unit operative to iteratively processing the output pattern by performing a convergent 
reverberation to yield a resonance pattern representative of resonances between a predetermined set 
of selected Quantum Expressor Functions and the output pattern until a predetermined degree of 
convergence is achieved between the resonances found in the resonance pattern and resonances 

25 expected for the set of mutations. 

26. The system of claim 23 further including a image preconditioning unit operative to 
perform one or more of the functions of: 

refocusing the output pattern to yield a re-focused output pattern; 
30 cross-correlating the re-focussed output pattern; 

applying a local maxima filter to the correlated re-focussed output pattern to yield a 
maximized output pattern; 

re-scaling the maximized output pattern to yield a uniformly re-scaled output pattern; and 
re-scaling the uniformly re-scaled output pattern to amplifying local boundaries therein to 
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27. The system of claim 23 wherein the resonance pattern interpretation unit is 
operative to map the confirmed mutations to known diseases associated with the pre-selected set of 

5 known mutations to identify diseases, if any, indicated by the biological sample. 

28. The system of step 23 further including a Quantum Expressor Function generator 
operative to: 

input a set of mutation signatures representative of the set of known mutations; 
1 0 input a representation of a biochip microarray pattern layout for the biochip from which the 

output pattern is generated; 

generate a set of resonant interaction parameters representative of mutation pattern 
interactions between elements of the biochip including interactions from a group including 
element-to-element interactions, element-to-ensemble interactions, ensemble-to-element 
15 interactions, and ensemble-to-ensemble interactions; and 

generate the set of nonlinear Quantum Expressor Functions from the set of resonant 
interaction patterns. 

29. The system of claim 23 wherein the output pattern is a quantized output pattern. 
20 30. The system of claim 29 wherein the quantized output pattern is a dot spectrogram. 

31. The system of claim 29 wherein the biological sample is selected from a group 
consisting of a DNA, RNA, protein, peptide-nucleic acid (PNA) and targeted nucleic amplification 
(TNA) samples. 

32. A system of generating a set of nonlinear Quantum Expressor Functions based upon 
25 a set of mutation signatures representative of a set of known mutations and based upon a biochip 

microarray pattern, said system comprising: 

a resonant interaction parameter generator operative to generate a set of resonant 
interaction parameters representative of mutation pattern interactions between elements of the 
biochip including interactions from a group including element-to-element interactions, element-to- 
30 ensemble interactions, ensemble-to-element interactions, and ensemble-to-ensemble interactions; 
and 

a Quantum Expressor Functions generator operative to generate a set of nonlinear Quantum 
Expressor Functions from the set of resonant interaction patterns. 
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